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Abstract. We have calculated constraints on the evolution of the equation of state 
of the dark energy, w(z), from a joint analysis of data from the cosmic microwave 
background, large scale structure and type-la supernovae. In order to probe the time- 
evolution of w we propose a new, simple paramctrization of w, which has the advantage 
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available. Furthermore it is well behaved in all asymptotic limits. Based on this 
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not have a better goodness-of-fit. Rank correlation tests on supernova data also do 
not show any need for a time- varying w. 
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1. Introduction 

The discovery in 1998 [1, 2] that the universal expansion is currently accelerating is one 
of the most spectacular result in cosmology from the past decade. The finding has since 
been confirmed by observations of the Cosmic Microwave Background (CMB) [3,4] and 
the large scale structure (LSS) of the universe [5-7] 

One possible explanation is that the energy density of the universe is dominated 
by dark energy with a negative equation of state. The simplest possibility is the 
cosmological constant which has P = wp with w — — 1 at all times. However, since 
the cosmological constant has a value completely different from theoretical expectations 
one is naturally led to consider other explanations for the dark energy. 

A light scalar field rolling in a very flat potential would for instance have a strongly 
negative equation of state, and would in the limit of a completely flat potential lead to 
w — — 1 [8-10]. Such models are generically known as quintessence models. The scalar 
field is usually assumed to be minimally coupled to matter, but very interesting effects 
can occur if this assumption is relaxed (see for instance [11]). 

In general such models would also require fine tuning in order to achieve fix ~ fi mi 
where fix and fl m are the dark energy and matter densities at present. However, by 
coupling quintessence to matter and radiation it is possible to achieve a tracking behavior 
of the scalar field so that fix ~ flm comes out naturally of the evolution equation for 
the scalar field [12-20]. 

Many other possibilities have been considered, like /c-essence, which is essentially 
a scalar field with a non-standard kinetic term [21-27]. It is also possible, although 
not without problems, to construct models which have w < — 1, the so-called phantom 
energy models [28-42]. 

Finally, there are even more exotic models where the cosmological acceleration is 
not provided by dark energy, but rather by a modification of the Friedman equation due 
to modifications of gravity on large scales [43,44]. 

Given this plethora of different possibilities and since we have no fundamental 
theory available for calculating it appears that w should be treated as an effective 
parameter only. In many models w is also changing with time, typically going from one 
asymptotic limit for z — > oo to another for 1 + z — > 0. 

The simplest parametrization is w = constant, for which constraints based on 
observational data have been calculated many times [45-48]. However, as the precision 
of observational data is increasing is it becoming feasible to search for time variation in 
w. 

In the present paper we propose a very simple parametrization for w which allows 
us to treat almost all models for dark energy currently on the market. Furthermore it 
is straightforward to extend our parametrization when better data becomes available, 
and the parametrization is well-behaved in all relevant asymptotic limits. The paper is 
structured as follows: In section 2 we describe our parametrization and its relation to 
supernova luminosity distances and CMB angular distances. In section 3 we discuss how 
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a numerical likelihood analysis has been performed with the most recent cosmological 
data, and in section 4 we discuss our results. We also discuss how to extend our current 
parametrization to provide a more refined description of w(z) when more data becomes 
available. Finally, section 5 contains a conclusion. 



2. Luminosity distance, angular distance, and the parametrization of w 



2.1. Supernova luminosity distances 

Type la supernova (SNIa) observations provide the currently most direct way to probe 
the dark energy at low to medium redshifts since the luminosity-distance relation is 
directly related to the expansion history of the universe. 
The luminosity distance di is given by 

' (l + z) 7 ^ m sm(^hl), Vl k <0 

d L = \(l + z)I, Q k = (1) 
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Here, w(z) is an arbitrary function of redshift. Putting w = constant = — 1, we have 
Qx = and the luminosity distance relation can be used to constrain the values of 
Q m and Q\. Using current Sla data, Q\ = can be ruled out at very high confidence 
level [49,50] (note that one marginalizes over an arbitrary multiplicative factor in dL 
and that the Hubble parameter H is thus treated as a free parameter). 

In principle, it should also be possible to constrain w(z) using SNIa data. In 
practice however, this is quite complicated since the luminosity distance depends on 
w through multiple integrals causing degeneracies between parameters [52]. Also, the 
sheer number of parameters to fit in relation to the quantity and quality of current 
data makes this difficult to do. As a simple rule of thumb, it is hard to constrain more 
than two parameters using current SNIa data. A common approach is to assume a flat 
universe, i.e., Qx = 1 — and fit Q m and w = constant. The flat universe assumption 
can be justified from a theoretical point of view as a consequence of inflation with 
observational support from CMB angular scale measurements. This approach gives an 
limit of w <^ — 0.5 using the "gold" sample of 157 SNIa compiled in [50]. Combining 
the SN data with CMB and LSS measurements yields w = — 1.08±[j;?8, consistent with 
a cosmological constant. 

The next step is to constrain the time variation in the equation of state parameter. 
One approach is to express w as a first order Taylor expansion in z (usually around 



Cosmological constraints on the dark energy equation of state and its evolution 



4 



^ = o)|| 

w(z) = w\ z=0 + — z. (2) 
dz z=o 

The advantage of this method is the simplicity; any value of dw/dz(z = 0) ^ 
would indicate a non-constant equation of state parameter. The obvious disadvantage 
is the divergence of w at high redshifts that makes it practically useful for SN data 
only^f. Assuming a flat universe and a prior of Q m = 0.27 ± 0.04 (derived assuming a 
cosmological constant), Riess et al. obtains w\ z=0 = — 1.31±q;|| and dw/dz(z = 0) = 
1.48±o;|q with the case of a cosmological constant within the (joint) 68 % confidence 
level. 

In additional to the constant and linear models of w there exists an ever increasing 
number of models for w that can be broadly classified into parametric models [53-61] 
and non-parametric models [63-67]. 

There are large differences in how to include information from CMB and LSS data, 
the choice of SN data sets, the marginalization over the SN magnitude zero point, the 
use of ejxternal priors in Q m and H etc. In spite of this, the consensus is that current 
data are consistent with dark energy in the form of a cosmological constant at the 2a 
level. The best fit model for a time varying equation of state tend to have lower values 
of w at low redshift. Whether this is a real effect or an artefact of the specific model 
employed in the fit is currently a subject of discussion [59, 68], but the fact remains that 
any indications of a possible time evolution in w are still statistically insignificant. 

In order to be able to constrain the temporal variation of w, it is crucial to combine 
different observational data in a self consistent way. For example, using priors from 
CMB and LSS derived assuming a cosmological constant when constraining w using 
SN data may lead inconsistencies. However, constraining w over a redshift range of 
< 2^1000 puts high demands on the parametrization (or non-parametric model) of 
w. 

2.2. Parametrization of w(a) 

In order to describe a possible time evolution of w(z), we propose the following very 
simple approximation 

w(a) = WqWx— — — , 3 

a q wi + alw 

where a is the scale factor normalized to the value of one today, i.e, a = (1 + z)~ l and 
wo,wi, a s , and q are constants, wo and w\ describe the asymptotic behaviour of w 

w(z) -> I (4) 
I Wi for 1 + z — > oo 

|| Of course, w can be expanded to higher order with the price of adding more parameters to the fit. 

A similar parametrization that avoids this behaviour at high redshift is w(a) = tu| a =i + w a (l — a) 
[69]. However, that parametrization diverges at a — > oo. 
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Figure 1. w(a/a s ) for different values of q. The values shown are q = 0.5, f , 2, 5, 10 
in order of increasing steepness. 

The two additional parameters a s and q describe the scale factor at changeover and 
the duration of the changeover in w respectively. In figure 1 we show the evolution of 
w as a function of redshift for different valus of q ranging from 0.5 to 10, the lower the 
value of q, the longer the duration of the changeover in w. In this plot, we have fixed 
Wo = —0.5 and ui\ = —1.5. 

Advantages of this parametrization include the simple and intuitive role of each 
parameter as well as the large flexibility that aids in avoiding problems such as "sweet 
spots" and fixed asymptotic values for w. 

2.3. CMB angular scales 

Changing the equation of state of the dark energy while keeping all other parameters of 
the model fixed has two distinct effects: a) It changes the magnitude of the IS W effect 
at low /, b) It introduces a linear shift in the angular position of the CMB features by 
a factor F, given by [48] 

F = Sl% 2 fitful + Zf + (1 - n m )(l + z )3(H-W)]-l/2. (5) 

Jo 

In figure 2 we show the value of F relative to that of an Q, m = 0.3 flat ACDM 
model. For each value of (wq,Wi), a s and q have been chosen so that F is as close to 
Fq = 1.61 as possible. From the figure it can be seen that a significant degeneracy in 
(wq,Wi) can be expected for CMB data. The reason simply is that F is an integral 
along a line of sight so that changing w by different amounts at different epochs can 
yield the same integrated result. CMB does not have access to the instantaneous value 
of w at each given redshift point, only to the integral F. 
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Figure 2. Contourplot of the function F. For each value of (wo,w\), a s and g have 
been chosen so that F is as close to F = 1.61 as possible. 

3. Numerical analysis 

Since it is not possible to obtain useful limits for all four parameters of our equation of 
state parametrization in equation 3 and we are primarily interested in exploring evidence 
for and/or a changing w, we employ the following approach: For fixed values 

of Wo and w\ we marginalize over the other cosmological parameters. As our framework 
we choose the minimum standard model with 6 parameters: Q m , the matter density, 
the curvature parameter, fl b , the baryon density, H , the Hubble parameter, and r, the 
optical depth to reionization. The normalization of both CMB and LSS spectra are 
taken to be free and unrelated parameters. The priors we use are given in Table 1. 

Table 1. Priors on cosmological parameters used in the likelihood analysis. 



Parameter Prior Distribution 



il = O m + tt x 1 Fixed 

h 0.72 ± 0.08 Gaussian [74] 

n b h 2 0.014-0.040 Top hat 

n s 0.6-1.4 Top hat 

t 0-1 Top hat 

Q Free 

b Free 



Likelihoods are calculated from x 2 so that for 1 parameter estimates, 68% confidence 
regions are determined by A% 2 = % 2 — Xo — 1) an d 95% region by A% 2 = 4. Xo is X 2 
for the best fit model found. In 2-dimensional plots the 68% and 95% regions are 
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formally defined by A% 2 = 2.30 and 6.17 respectively. Note that this means that the 
68% and 95% contours are not necessarily equivalent to the same confidence level for 
single parameter estimates. 

3.1. Supernova luminosity distances 

We perform our likelihood analysis using the "gold" dataset compiled and described in 
Riess et al [50] consisting of 157 SNIae + using a modified version of the SNOC package 
[51]. 

3.2. Large Scale Structure (LSS). 

At present there are two large galaxy surveys of comparable size, the Sloan Digital Sky 
Survey (SDSS) [6, 7] and the 2dFGRS (2 degree Field Galaxy Redshift Survey) [5]. Once 
the SDSS is completed in 2005 it will be significantly larger and more accurate than the 
2dFGRS. In the present analysis we use data from SDSS, but the results would be almost 
identical had we used 2dF data instead. In the data analysis we use only data points 
on scales larger than k = 0.15/i/Mpc in order to avoid problems with non-linearity. 

3.3. Cosmic Microwave Background. 

The CMB temperature fluctuations are conveniently described in terms of the spherical 
harmonics power spectrum Cf T = (|a/ m | 2 ), where ^-(#,0) = Y,im a imYim{Q , 0) • Since 
Thomson scattering polarizes light, there are also power spectra coming from the 
polarization. The polarization can be divided into a curl-free (E) and a curl (B) 
component, yielding four independent power spectra: Cf T , C EE , Cf B , and the T-E 
cross-correlation Cf E . 

The WMAP experiment has reported data only on Cf T and Cf E as described in 
Refs. [3,4,70-72]. We have performed our likelihood analysis using the prescription 
given by the WMAP collaboration [4, 70-72] which includes the correlation between 
different Cj's. Foreground contamination has already been subtracted from their 
published data. 

4. Results 

In Figs. 3-5 we show results of the likelihood analysis. For each (w ,Wi) we have 
marginalized over all other parameters. 

From the figure including CMB and LSS data only it is clear that there is an almost 
complete parameter degeneracy with the same shape as predicted in Fig. 2, i.e. it follows 
the degeneracy inherent in the angular distance integral. With the present data there is 

+ Note that the electronic table in ApJ is one SN short. To get the full data set, use the table in 
astro-ph version. The missing SN is SN1991ag. 
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Figure 3. 68% (dark) and 95% (light) likelihood contours for WMAP and SDSS data 
only. 




Figure 4. 68% (dark) and 95% (light) likelihood contours for SNI-a data only. 

no way of discriminating between models with strong time variation of w, and models 
with constant w. 

When only SN data is used there is also an almost perfect degeneracy which comes 
from the integral I in Eq. (1). Again, there is no means of discriminating between a 
time varying and a constant w. 

However, once SN, CMB and LSS data are combined this degeneracy is partially 
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Figure 5. 68% (dark) and 95% (light) likelihood contours for WMAP, SDSS and 
SNI-a data. 

broken. From Fig. 5 it is clear that the best fit has w± > wq, meaning that w is decreasing 
with time. The reason that Fig. 5 is not a superposition of Figs. 3 and 4 is that the best 
fits for the individual subsets of data require different values of Q m , a s and q. In Table 2 
we tabulate x 2 f° r various best fit models, compared with the corresponding best fits 
for ACDM models (wq = W\ = —1). The best fit model to the combined CMB, LSS and 
SNI-a data has \ 2 — 1623.1, which should be compared to the 1512 degrees of freedom 
in the fit. The goodness-of-fit (GoF)* of this model is 0.0236, indicating that only 2.36% 
of all randomly generated data sets based on the same underlying cosmological model 
would produce a worse fit than the one observed. The main cause of this bad fit is the 
WMAP data, a fact which is well known, and which is presumably due to incomplete 
foreground removal. 

At first sight the models with varying w seem to better fit the data, disfavouring a 
pure cosmological constant. However, four new parameters are added to the model with 
time varying w. The GoF for the best fit ACDM model is actually 0.0239 {\ 2 = 1626.9 
for 1516 degrees of freedom), which is better than for the time- varying w model. 
Therefore there is no compelling evidence for a time-variation of w in the present data. 
We elaborate more on this point in section 4.4 where we perform non-parametric rank 
correlation tests on the SN data. 

Our best fit model has parameters Wi = —0.4, w = —1.8, Q m = 0.38, a s = 0.50, 
and q = 3.41. In figure 6 we show w(z) for this model. Evidently the data favour a 
strongly time-dependent w at present. We note that this result is similar to that of [53]. 

* Defined as GoF = T(i'/2, \ 2 /2) /T(i//2), where T(zy/2,x 2 /2) is the incomplete gamma function and 
v is the number of degrees of freedom. 
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Figure 6. Evolution of w in the best-fit model found. 

In this work w is also parameterized as evolving from one asymptotic limit to another, 
albeit with a different parametrization. Their best fit model is found to have w± = —0.7, 
wo = —2.0, using slightly different observational data than is used in the present work. 

Table 2. x 2 for the best fit model, compared with ACDM, for different data sets 



Data used 


best fit 


ACDM 


CMB+LSS 


1446.7 


1447.3 


SNI-a 


173.8 


177.1 


CMB+LSS+SNI-a 


1623.1 


1626.9 



4-1. A constant equation of state 

The simplest possibility for dark energy is that w is constant. Constraints on this type 
of model can be found by setting w — w±, effectively taking line sections through the 
(w ,Wi) likelihood plots. The resulting ID x 2 curves can be seen in figure 7. From the 
combination of all available data we find that —1.34 < w < —0.79, in good agreement 
with other recent analyses [49,50]. 

4-2. Parametrization in terms ofw(z = 0) and dw/dz(z = 0) 

In fact our parametrization can easily be reduced to a model where the fit is done for 
w(z = 0) = u>* and dw/dz(z — 0) = w'^. 
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Figure 7. Likelihood functions for the case of wq — w%. The full (red) line is for 
WMAP+SDSS, and the dashed (blue) for WMAP+SDSS+SNI-a data. 



To do a one dimensional fit over w#, one simply rewrites W\ as 
w*a q s w 

W\ = ; (6) 

w (l + at) - iu* 

When the likelihood analysis is done, wq is a free parameter and W\ given by the 
relationship above. The result of this exercise is shown in Fig. 8. From this it can 
be seen that the best fit is at u>* = —1.43 as described above, and that io* = — 1 is 
disfavoured at about 2a, which is also in accordance with our previous finding. The la 
allowed interval is u>* = — 1.43loil- While the best fit value is very close to what is 
obtained by Alam et al. [58,59], and also by Riess et a!., the error bars are different 
from those of Riess et al. The reason is that we add other observational data which 
shift the best fit value, and also that differences between our parametrization and the 
Taylor series used in other analyses increase when moving away from the best fit values. 
Basically our parametrization has more freedom to fit data because there are additional 
free parameters. This in turn means that the one-parameter error bars increase. 

Equally, it is possible to do a parametrization in terms of w 1 ^ instead in the following 
way. First, the derivative of w with respect to a is 



, dw 
da 



_ w w 1 qa q a (wi - w ) 
=0 (wi + aiw ) 2 

We now introduce a new variable x = a q s . By doing this, we can solve the above equation 
for q 

, Qi + XWg) 2 

V = w * 7 y 8 

W VJiX{W — Wi) 
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Figure 8. \ 2 as a function of the present day value of w for the WMAP+SDSS+SNI-a 
data. 

In the likelihood analysis, the free parameters are then wq, ui\ and x, whereas q is given 
by the above relation. The result of this analysis is shown in Fig. 9 

Formally, the la allowed interval is w'^ = 1.0iS;g. However, as is evident from the 
figure the likelihood is far from Gaussian, and no useful constraints can be placed at 
the 2<T level. 



4-3. Comparison with analyses based on Taylor series expansions 

Tentative evidence for evolution of w has been found in several other analyses of 
supernova data. Most analyses of SN data use a Taylor series expansion of w around 
z = because no high redshift data is used. One recent example is the work of Alam 
et al. [58, 59] where the parametrization 

Pde = A + A 1 (l + z) + A 2 (l + zf (9) 

is used. This can be translated into a relation for w 

1 + z A 1 + 2A 2 (l + z) 
'"'W— 1 + — 4, + + ,) + *(! + »)' ' <10) 
where A = 1 — Q m — A\ — A 2 . 

The best fit found in that work is fl m = 0.3, Ai = —4.36, A 2 = 1.83 (slightly better 
fits are found for higher Q m , but the corresponding Ai and A 2 are not tabulated). This 
leads to 

dw 



w\ z=0 = -1.33 — 
dz 



= 1.07 (11) 

z=0 
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Figure 9. \ 2 as a function of the present day value of dw/dz for the combined 
WMAP+SDSS+SNI-a data. 



Our results cannot be compared on global scale because the Taylor series expansion 
quickly breaks down (in fact it cannot generally be expected to be valid even at z ~ 1). 
The only direct comparison available is in terms of the present w and its derivative. 
As described above, we find that the best fit model has 

dw 



w\ z=0 = -1.4 

dz 



= 1.0 (12) 

=o 

which is almost identical to what is found by Alam et al. 

As discussed previously the reason that dw/dz ^ is a better fit than a constant 
w is exactly the supernova data, and therefore our results can be expected to be very 
similar in the one limit where they can be directly compared. 



4-4- SNIa data rank correlation test 

Since the SNIa data are best fit by a changing w and there are different opinions about 
the significance of the effect and the dependence on the specific dark energy model, 
we have also investigated the need for additional parameters in the luminosity distance 
relation by studying the correlation between the observed magnitude and the redshift 
after subtracting the best fit cosmology magnitudes. This can be done in an non- 
parametric manner using the Spearman or Kendall rank correlation tests. The two- 
sided significance level, p (ranging from to 1), of the deviation from zero correlation 
indicates whether the model is capable of representing the complexity of the actual data 
or not. A small value of p indicates a significant correlation or anticorrelationjj. For an 

H Note that the zero point magnitude is irrelevant for both the likelihood analysis and the correlation 
tests. 
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Einstein-de Sitter universe, we get p<^ 10~ 19 , indicating - as expected - that this simple 
model is not able to mimic the observed magnitudes over a wide redshift range. Fitting 
Q m for a zero cosmological constant universe, we get p ~ 0.17 — 0.19 (for Q m = 0), 
showing a better but still unsatisfactory fit. Including a cosmological constant with or 
without a flat universe assumption removes any evidence of correlation (see table 3) 
and disfavouring the need for more complexity in the dark energy model. We have also 
simulated a sample 100 SNIa data sets with the same number and redshift distribution 
of SNe as the Riess gold sample assuming a cosmological constant flat universe with 
VL m = 0.3 and a dispersion in SN magnitudes of 0.3 mag. Subtracting the best fit 
cosmological constant flat universe and performing the rank correlation tests on each 
simulation yields 

Vt m = 0.31 ±0.05 
p s =0.47 ±0.29 

p K = 0.47 ±0.30. (13) 

Comparing with the results in table 3 it is clear that the real data actually have a smaller 
deviation from zero correlation than is expected for a cosmological constant universe, 
thus weakening the case for a evolving w. 

This result is completely consistent with the fact that the Goodness-of-Fit for the 
combined CMB, LSS, and SNI-a analysis is slightly worse for the time- varying w models 
than for pure ACDM models. 

Table 3. Result from likelihood analysis and correlation test, ps and px are the 
two-sided significance level (ranging from to 1) of the deviation from zero for the 
Spearman and Kendall rank correlation tests. A small value of p indicates a significant 
correlation or anticorrclation. 



Cosmology 


Best fit parameters 


Ps 


PK 


Einstein-de Sitter 








io- 19 





Matter dominated 




= 




0.20 


0.17 


fl m + f?A = 1 




= 0.31 




0.88 


0.86 


VL m and f^A 




= 0.46, 0\ = 


0.98 


0.88 


0.86 


Flat with const, w 




= 0.49, w = - 


-2.4 


0.88 


0.89 


Flat with linear w 




= 0.48, to, = 


-2.5, < = 2.0 


0.98 


0.98 



4-5. Extending the parametrization 

The parametrization used in the present analysis already involves more parameters than 
can be unambiguously fitted with present data. Nevertheless it it still interesting to see 
how the parametrization naturally extends to a larger number of parameters. 

If we still use the assumption that the transition is between two asymptotic values 
of w, the next natural parameter to add is a skewness in the transition profile. This can 
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Figure 10. w(a/a s ) for different values of r. Other parameters have been fixed to 
wo = —1.5, w\ = —0.5, q = 3. The full line is for r = 0, the dashed for r = 2 and the 
dotted for r = —2. 



for instance be incorporated with the parametrization 

a q+(a/a a ) r -l , a q 

w(a) = w Wi — -j—. 1 — , (14) 

where r is a new parameter. For negative r the transition is more rapid at early times, 
whereas for positive r it is more rapid at late times. In fig. 10 we show the effect of 
adding r as a new parameter. 

Another possibility is that there are several transitions in w during the evolution 
of the universe. This is also simple to incorporate. The initial transition is described by 

Wi(a) =w w 1 a +CLs — . (15) 
a q uii + aiwo 

To add a second transition, simply take 

w(a) = Wi(a)w 2 — r s2 . r , (16) 

a r w 2 + a r s2 Wi(a) 

where a s2 is the scale factor at the transition w — > wi, and r describes how rapid the 
transition is. In figure 11 we show a model which has the second transition at a s2 = 20a s , 
q = 5, r = 2, w\ = —0.5, w = —1.5, and w 2 = —1.7. By this iterative procedure, more 
and more steps can be added to the parametrization. 

Finally, we note that our parametrization is well-behaved only if w < 0, even 
though w can be arbitrarily close to 0. If there is a transition from w > to w < 
there is a singularity in w. However, this can easily be cured by taking instead 
w(a) = (a q w + ajwi)/(a q + af). 
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Figure 11. w{a/a s ) for a model with two transitions. 
5. Conclusion 

We have analysed present cosmological data, using an extension of the current minimal 
cosmological model to allow for time-variation in the dark energy equation of state, w. 

The result of the likelihood analysis is that the data favour a w which is currently 
undergoing a rapid transition towards a more negative value. In our parametrization 
where w varies between two asymptotic values, the best fit is w{z — > oo) = —0.4, 
and w(l + z — > 0) = —1.8. This result agrees well with other analyses using different 
parametrizations. 

The la allowed interval for the present value of w is w(z = 0) = — 1.43zto!38 an d for 
its derivative dw/dz(z = 0) = l.Ol^g. 

Even though pure ACDM models have higher y 2 values, they also have more fitting 
parameters. When the Goodness-of-Fit is calculated for both types of models, the best 
fit ACDM model has GoF = 0.0239, whereas the time-varying w model has GoF = 
0.0236. This indicates that there is no real evidence for a time variation of w in the 
present data. This finding is corroborated by rank correlation tests on the SNI-a data. 
From Spearman and Kendalls tests we again find that time-varying w models do not 
produce a statistically significant improvement in the fit to data. 

Finally we have discussed in some detail how to extend our parametrization of w 
to allow for a more refined description once better data becomes available. 
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